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Abstract. We consider limiting Carleman weights for Dirac operators 
and prove corresponding Carleman estimates. In particular, we show 
that limiting Carleman weights for the Laplacian also serve as limiting 
weights for Dirac operators. As an application we consider the inverse 
problem of recovering a Lipschitz continuous magnetic field and electric 
potential from boundary measurements for the Pauli Dirac operator. 



1. Introduction 

In this article we consider Carleman estimates for Dirac operators, with 
applications to inverse problems. There is an extensive literature on Car- 
leman estimates and their use in unique continuation problems (see for in- 
stance tlSj). However, Carleman estimates have also become an increasingly 
useful tool in inverse problems, we refer to [H] for some developments. 

The recent work |17] shows the importance of Carleman estimates in the 
construction of complex geometrical optics (CGO) solutions to Schrodinger 
equations. CGO solutions have been central in establishing unique identifia- 
bility results for inverse problems for elliptic equations (starting with [3] and 
[29] , see surveys [30] , [H] ) . The approach of [TT] provides a general method 
of constructing CGO solutions, and has led to new results in inverse prob- 
lems with partial data [7], [18], determination of inclusions and obstacles 
[33] . [13] ■ [32], [28], and also anisotropic inverse problems [6j. All the re- 
sults listed are concerned with, or can be reduced to, Schrodinger equations 
whose principal part is the Laplacian. 

Our aim is to extend the Carleman estimate approach of [T7] to Dirac 
operators and their perturbations. By a Dirac operator we mean a self- 
adjoint, homogeneous, first order N x N matrix operator Pq with constant 
coefficients in R", whose square is —A. Examples include Dirac operators 
on differential forms in R", and the Pauli Dirac in R^ given by the 4x4 
matrix 
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where D = —iSJ and a = (o"i, o"2, (T3) is a vector of Pauh matrices with 




The first step is to prove Carleman estimates for Dirac with a special 
class of weights, called limiting Carleman weights. These weights were in- 
troduced in [17] for the Laplacian (see Section [2] for the definition) , and in 
this case include linear and logarithmic functions. See [6] for a complete 
characterization. One may think of a limiting weight as a function ip such 
that the Carleman estimate is valid both for 99 and —ip. A typical result for 
Dirac is as follows. Here and below || • || = || • and A< B means that 

A < CB where C > is a constant independent of h and u. 

Theorem 1.1. Let Q C R", n > 2, be a bounded open set, let Pq be a Dirac 
operator as above, and let V G L°°{^)^^^ . Let cp he a smooth real-valued 
function with Vip 7^ near 0,, and assume that ip is a limiting Carleman 
weight for the Laplacian. If /i > is sufficiently small, then one has the 
estimate 

h\\u\\ + h^\\Vu\\ < We'^^^iPoihD) + hV)e-'^/^u\\, 
for any u G C^{n)^ . 

We give two different proofs of such estimates. The first proof uses the fact 
that Dirac squared is the Laplacian, and gives the estimate for any limiting 
Carleman weight for the Laplacian. This proof, however, is based on symbol 
calculus and does not give Carleman estimates with boundary terms, which 
were crucial in the partial data result of [17]. We will present another proof, 
which uses just integration by parts and gives a simple Carleman estimate 
with boundary terms. The second proof is valid for the linear weight. 

The next step is to construct CGO solutions to Dirac equations by us- 
ing the Carleman estimate. These solutions are then used to prove unique 
identifiability of the coefficients from boundary measurements. We carry 
out this program for an inverse problem for the Pauli Dirac operator with 
magnetic and electric potentials, which we set out to define. 

Let Pq be as in ([T]) and consider the operator 

Cv = PoiD) + V, (2) 
where the potential V has the form 

V = Po{A)+Q, (3) 

with Q = ("V^ J, J, A = (ai, 02,03) G t^i'~(l^;R3), and q± G W^^°-{n). 

We consider Cy acting on 4-vectors u = (u- ) where u± G L^(ri)^. It is 
shown in [22j that the boundary value problem 

J Cyu =0 in Q, 
\ u+ = f on d^}, 

is well posed if has smooth boundary and is not in the spectrum of 
Cy There is a unique solution u G Ti.{0,)'^ for any / G /i(5ri), where 
n{n) = {u£ L2(0)2 ; a • L>m G L2(17)2} and h{dn) = nln)\an- We refer to 
[22j for more details and a description of the trace space h{di}). 
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There is a well-defined Dirichlet-to-Dirichlet map 

Kv : h{dQ) ^ h{dQ), f ^ u^lan. (4) 

Thus, Ay takes ti+ to u_ on the boundary, where u is the solution of the 
Dirichlet problem above. This map is invariant under gauge transformations: 
if A is replaced by ^ + Vp where p\qq = 0, the map Ay stays unchanged 
and so does the magnetic field V x ^. 

More generally, if Q is any bounded open set in R", we define the Cauchy 
data set 

Cy = {{u^\qq, u-lgn) ; u G H{^1)'^ and £yu = in ^l}. 

Here the traces are taken in the abstract sense as elements of Tl{i^)/ Hq{Q). 
If is smooth and is outside the spectrum, Cy is just the graph of Ay. 
Also here Cy is unchanged in the gauge transformation A ^ A + Vp if 
plan = 0. 

We consider Cy as the boundary measurements, and the inverse prob- 
lem is to determine the magnetic field V x yl and the electric potentials 
q± from Cy. This inverse problem, and the correspoding inverse scattering 
problem at fixed energy, have been studied in pjj, [12j, [22], [20] with vary- 
ing assumptions on the potential. The next uniqueness result was proved 
for smooth coefficients in [22] by reducing to a second order equation and 
using pseudo differential conjugation. We establish this result for Lipschitz 
continuous coefficients by working with the Dirac system directly using the 
Carleman estimate approach. In particular, the proof avoids the pseudodif- 
ferential conjugation argument. 

Theorem 1.2. Let C be a bounded C^ domain, let Aj e W'^'°°{n; R^), 
and let q±j G W^''^{n), j = 1,2. If Cy, = Cy^, then V x = V x ^2 and 
q±,i = Q±,2 in 0. 

The first step in the proof is a boundary determination result, which again 
was proved in [22] for smooth coefficients and domains by pseudodifferential 
methods but which is needed here in the nonsmooth case. We use argu- 
ments of [2], [3] given for scalar equations, and construct solutions to the 
Dirac equation which concentrate near a boundary point and oscillate at the 
boundary. These solutions can be used to find the values of the tangential 
component ^tan = A — {A ■ u)u of A and of q± at the boundary. Here v is 
the outer unit normal to d^. 

Theorem 1.3. Let C R^ be a bounded C^ domain and A, q± £ W^'°°{n). 
Then Cy uniquely determines ^tanlan and q±\dn- 

Besides their independent interest, a major motivation for studying Dirac 
operators comes from inverse problems for other elliptic systems. These 
problems have often been treated by reducing the system to a second order 
Schrodinger equation, such as with the Maxwell equations [25j or linear 
elasticity [23], [21], [TT]. These reductions are however not without problems, 
and for instance counterparts of the partial data result in [17] for systems 
may be difficult to prove in this way. For the Maxwell equations, there is 
another useful reduction to a 8 x 8 Dirac system (see [25] ) which may be 
more amenable to a partial data result. 
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Another benefit of working with first order systems directly is that the 
reduction to a second order system may require extra derivatives of the co- 
efficients. For instance, the reduction of Maxwell equations to a Schrodinger 
equation in |25j requires two derivatives of the coefficients, while the Dirac 
system only requires one derivative. Also, the Carleman estimates are valid 
for low regularity coefficients. Theorem II. 2t whose proof uses the first order 
structure and Carleman estimates, may be thought of as an example result 
with improved regularity assumptions. 

For other work on Carleman estimates and unique continuation for Dirac 
operators, we refer to [l], [16], [21]. Also, while writing this article, we 
became aware of the work of M. Eller [9], which also considers Carleman 
estimates for first order systems via integration by parts arguments. 

The structure of the article is as follows. In Section 2 we give different 
proofs of Carleman estimates for Dirac with limiting weights. Section 3 gives 
a construction of CGO solutions for ([2j) with smooth coefficients, and the 
case of Lipschitz coefficients is considered in Section 4. The last two sections 
prove the uniqueness results for the inverse problem. Theorems 11.21 and 11.31 

2. Carleman estimates 

Let Pq{^) he an N X N matrix depending on ^ G R", such that each 
element of fb(0 is of the form a ■ ^ where a G C". We assume the two 
conditions 

Poio' = eiN, Poio* = Poio, 

for all 1^ G R". These conditions imply that the operator PQ^hD), where 
/i > is a small parameter, is a self-adjoint semiclassical Dirac operator: 

PoihDf = {-h^A)lN, PoihD)* = Po{hD). 

We may as well define Po(0 for ^ G C". Then Po(0* = PoiO, but one still 
has PoiO'^ — C'^^N- This and the linearity of Pq imply 

PoiC)Po{0 + PoiC)Po{C) = 2iC-OlN forCeeC". (5) 

In this section 0, will be a bounded open set in R", where n > 2. We 
wish to prove a Carleman estimate, which is a bound from below for the 
conjugated operator 

Po,^ = e^/''Po{hD)e-^/\ 

It is well known that such bounds exist if the weight ip enjoys some pseu- 
doconvexity properties. On the other hand, in the applications to inverse 
problems one needs this estimate both for ip and — and the weight can 
only satisfy a degenerate pseudoconvexity assumption. This is the case for 
the limiting Carleman weights introduced in [T7], where also a convexifica- 
tion procedure was given for obtaining Carleman estimates for these special 
weights. 

We recall that <p £ C°°(i7; R) is a limiting Carleman weight for the 
Laplacian in the open set fi, where Q CC 0,, if V(/7 7^ and {a, 6} = 
when = 6 = 0, where a and b are the real and imaginary parts, re- 
spectively, of the semiclassical Weyl symbol of the conjugated scalar oper- 
ator 

g.p/h(_/^2^)g-vp/h^ Examples are the linear weight ip{x) = a ■ x where 
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a E R"", |a| = 1, and the logarithmic weight '^{x) = log \x — xq\ where 

One can deduce a Carleman estimate for the perturbed Dirac operator 
Po{hD) + hV directly from a corresponding estimate for the Laplacian. We 
use the semiclassical Sobolev spaces 



L2 



where (0 = (1 + I^P)^/^- 

The estimate for Laplacian is as follows. This was proved in |17j with a 
gain of one derivative, i.e. ||?x||j:^s+i controlled by the right hand side. We will 

scl 

give a slightly different proof using the Carding inequality to obtain a gain of 
two derivatives, which will give an improved estimate for Dirac. Below, we 
will use the conventions of semiclassical calculus. In particular we consider 
the usual semiclassical symbol classes 5"*, and relate to symbols a G 5™ 
operators A = Op f^{a) via Weyl quantization. See [5] for more details. Also, 
we write {u\v) = J^u-vdx and ||ii|| = (ulu)^^"^. 

Lemma 2.1. Let ip he a limiting Carleman weight for the Laplacian, and 
let = (/3 + be a convexified weight. Then for /i ^ e ^ 1 and s G R, 

-^||n|U.+2 < C\\e'^-/''(-h^A)lNe-^-/''u\\H^ (6) 

for ah u G C^{n)^. 

Proof. First extend if smoothly outside 0, so that ip = 1 outside a large 
ball. In this way one can use Weyl symbol calculus in R", and the choice of 
extension does not affect the final estimates since the differential operators 
are only applied to functions supported in Q. 

Let = e'^/''{-h'^A)e-'f/^, and write P^ = A + iB with A and B self- 
adjoint, so A = -h'^A - |Vy?p, B = VipohD + hDoVip. live C^{^) we 
integrate by parts to get 

\\P^vf = {{A + iB)v\{A + iB)v) = \\Avf + \\Bvf + {i[A, B]v\v). (7) 

Let the semiclassical Weyl symbols of A and B be given by a = — |V(/9p 
and b = 2'SJip-^. We recall the composition formula for semiclassical symbols 
p and q: the operator R = PQ has symbol 



Jj |a+/3| ('_-| \ |a| u 

The commutator is given by 

i\A,B\ = Ov^(}i{aM)- 

We now replace iphy ip^ = f{<p) where /(A) = A + j^- Let A,, and B^ 
be the corresponding operators with symbols and 5^, and let rj = f'{^p)(,. 
A computation in [T7] implies that whenever as{x,rj) = bs{x,rj) = one has 

{a„ rj) = 4n^)f'i^f\V^\^ + f'{^f{a, b}{x, 0- 
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For the following facts on symbols we assume that x is near Vt. The limiting 
Carleman condition and some algebra (see [T7], [18]) show that 

{ae,he}{x,r]) = Af"{Lp)f'{ipf\Vip\'^ + me{x)ae{x,ri) + ls{x,r])be{x,r]), 

where 



The symbol of i[A^, B^] is then given by 

h{ae, be} = /'((/?)^|Vv?|^ + hrueae + hlebe- 



The first term is always positive near il. since Vip is nonvanishing. To obtain 
a ifg^j bound we introduce the symbol 

de = Anipf\Vip\^+4+bi 

and we write the earlier identity as 

h{ae, be} = —de + hnieae + hlebe - —{al + bl). (8) 

Suppose that /i <C e <C 1. Since is elliptic and of order 2, there is a 
constant cq > 0, independent of e, such that 

de{x,r]) > co{ri)^, x near J7, 77 G R". 

The easy Garding inequality implies that for h small, 



{DeV\v)>^\\v\\l, , veC^m. 

scl 



On the other hand, we have the quantizations 
Op;,(a2) =Al^h^qx(x), 
Opi,ib^e) = B^ + h^Q2ix), 

Op^ (m^ae) = ^me o Ae + ^Ae ome + h^qsix), 
Ophihbe) = ^LeBe + ^BeLe + h^q^ix), 

where the qj are smooth functions which, together with their derivatives, 
are bounded uniformly with respect to £ near fl. It follows from ([S]) that 

h 

{i[Ae, Be]v\v) = — {DeV\v) + - [{AeV\meV) + {meV\AeV) 

+ {BeV\LeV) + {LeV\BeV)] - —{\\Aevf + \\Bevf ) - h^{q5v\v) 
> ? — 11^11^2 - Ch\\AeV\\ \\v\\ - Ch\\BeV\\ \\v\\jji 

Z £ scl scl 

_l!^(^\\A,vf + \\Bevf)-Ch^v\\\ 
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Recall that /i <C e <C 1, where e is fixed (depending on the constants C). 
Using the inequality \a(5\ < (5|ap + we obtain 

Now ([7]) shows that 

l|F...P>t^lMli.V (9) 

Finally, we need to shift ([9]) to a different Sobolev index to prove dH) . Let 
O CC with (f a limiting Carleman weight in 0,, and let x ^ C^(^) with 
X = 1 near Q. If u G C^{Q), we have the pseudolocal estimate 

Wil -x){hDyu\\H^^ <CMh^^\\u\\^, , for any a,/3,t,M. 



,cl 



This and ([9]) imply 

h\\{hDy+^u\\ = h\\{hDyu\\H2 

scl 

<h\\x{hDru\\H2+h\\{l-x){hDru\\H2^ 

scl scl 

<v^\\p^^ix{hDru)\\+h^\\{hDr+M- 

By the pseudolocal property we have 

\\[P^^,x]{hDru\\ < CMh^'\\{hDr^\\\, for any M, 

and since x[P(fie^ i^D)^] = hOpf^^S'^), we have for h <^ e <^ 1 that 

\\x[P^^,{hDr]u\\<h\\{hDru\\. 

One can check that all constants here can be chosen independent of e. From 
these estimates we obtain (j6]). □ 

Our first estimate for Dirac follows immediately. 

Lemma 2.2. Let V G L°°(r2)^^^ and let ip he a limiting Carleman weight 
for the Laplacian. If — 1 < s < 0, then for h small one has 

for all u G C^{n)^. 

Proof. Consider first Pq,^, = e'^'/^Po{hD)e-'^'/^. Now dSD implies 

A||^|| < C\\P^ u\\ , = C\\{hD)-'Po^^^Po,^M\Hi^, < C\\Po,^M\hi^, 

A / ^ scl scl 

(10) 

since {hD) ^Po,^ps is of order 0. We add the perturbation hV and use the 
estimate ||yn||H= < ll^^ll ^ ll""!! ^ ll^ll h-'^+i to get 

scl 

-^II-uIIh-s+i < C'(||(Po,¥>e + hV)u\\H^ + h\\V\\L°°\\u\\rTS+i). 

a/ £ scl ^"^^ scl 

Choosing e small enough we may absorb the last term to the left hand 
side. Since Po,^,+hV = e^^/^^e^/^{PQ{hD) + hV)e-^/^e-'^''/'^^ , we have the 
desired estimate. □ 

By using duality and the Hahn-Banach theorem in a standard way, one 
can convert this estimate (the case s = — 1) into a solvability result. 
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Proposition 2.3. Let V G L°°(0)^^^ and let (/J be a limiting Carleman 
weight for the Laplacian. If h is small, then for any / G L^(ri)^ there is a 
solution u G H^{Q)^ of the equation 



which satisfies < ||/||. 

scl 

In [T7] the partial data results for inverse problems were based on Carle- 
man estimates with boundary terms, that is, estimates proved for functions 
whose support may extend up to the boundary. The above proof of the Car- 
leman estimate for Dirac involved symbol calculus and estimates shifted to 
a different Sobolev index, which are problematic when applied to functions 
which are not compactly supported. 

We will present another proof of a Carleman estimate for Dirac, which 
involves only integration by parts and gives an estimate with boundary 
terms. The estimate is weaker than the earlier ones, and works for y = 0. 
Following [T7], we try to prove the estimate with boundary terms by writing 
the conjug ated Dirac Po,</p = e'P/^PQ{hD)e-^/^ as A + iB where A and B 
are self-adjoint, and by computing 

||Po,^nf = {{A + iB)u\{A + iB)u) 

= \\Au\\^ + \\Bu\\^ + {i[A^ B]u\u) + boundary terms. 

In |17j . which considered scalar operators, it was essential that the principal 
symbol i{ah — ha) of i[A,B] vanishes. Here however a = -Po(C) & = 
Po(V(/3) are matrices, and the principal symbol does not vanish. A part of 
it does vanish however: \iv = Vip/\Vip\ and a = ay + a^, where ay = Po(^||) 
with = • v)v, then 



This idea motivates the following proof. Here let {u\v)qq^ = Jq^ u ■ v dS and 
dn± = {xedn; ±Vip{x) ■ u{x) > 0}. 

Proposition 2.4. Let ^{x) = a ■ x where a G R", \a\ = 1. Let h be 
sufficiently small. Then for any u G C°°(r2)^, one has the Carleman estimate 

- hiid^^)u\u)dn. < \\e^^^PoihD)e-'p/''uf + hiid^ip)u\u)9n+ 

- ih{A\\u\Poiiy±)u)dn - ih{A^u\Po{iy\\)u)9Q, 

where A^^ = PQ{a)a ■ hD, z^y = (a • z^)a, and A± = PQ{hD) — Ay , i'± = — 

Proof. Let first cp G C~(J7;R) and V99 / on H. Let Pq,^ = A + iB he 
as above, so that A = PqQiD) and B = Pq{S/(p). Let v = V(p/\Vip\ = 
. . . , Vn)- We decompose A as A = A\\+ A±, where 



e'P/''{Po{hD) + hV)e-'P/^u = f in 17 



(11) 



aub - bau = (e • V^){Po{v)^ - Po{vf) = 0. 



n 



n 



A\\=Y,Poiv)vjhDj, A 



■± 



PoihD)-Y,Poiv)vjhDj. 



It holds that 



Po,^uf = {{A + iB)u\{A + iB)u) = \\A\\uf + \\{A± + iB)uf 

+ {A\\u\A±u) + {A±u\A\\u) + i{Bu\A\\u) - i{A\\u\Bu). 
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Use Po{v)'^ = I and integrate by parts to obtain 
i{Bu\A^^u) — i{A^^u\Bu) 

n 

= ^[i{PQ{VLp)u\Po{v)v.jhDju) - i{Po{v)vjhDju\Po{V^)u)] 

n 

= ^[^(^Iv'aij-^-Djn) - i{Lpx^hDju\u)] 

= -h{{dy(p)u\u)QQ, + /i((A(^)n|ti). 
For the terms involving A\\ and j4^, one has 
{A\\u\A±_u) + {A^u\A\\u) 

= {{A\A\\ + A\Ai_)u\u) + ih{A\\u\PQ{ui_)u)Qn + ih{Ax_u\PQ{u\\)u)QQ, 

where v\\ = {v ■ v)v and = — z^y . So far we have not used any special 
properties of (p. 

Now assume that (p{x) = q • x is the linear weight and choose coordinates 
so that f{x) = xi. Then t; = ei is a constant vector, and 

n 

^11 = Po{ei)hDi, A^ = Y,Po{ej)hDj. 

i=2 

We have ^| = A\\ and A\ = A^, and 

n 

A^A\\ + A\\A^ = Y,{Po{ej)Po{ei) + Po{ei)Po{ej))hDihDj. 

This vanishes by ([5]). 

From the above, we have for u G C°°(f2)^ the Carleman estimate 

\\Po,^uf = \\A\\uf + \\{A^ + iB)uf 

- h{{duip)u\u)aQ + ih{A\\u\PQ{h'i_)u)dQ, + ih{A^u\PQ{u\\)u)dn. 
The desired estimate follows. □ 

Remark. The above proof can be carried out for the convexified linear 
weight ip^ = if + where ip{x) = a ■ x, since u is a constant vector also 
in that case. In this way one can include a L°° potential in the Carleman 
estimate for Lp^, as well as the term on the left. However, due to the 

boundary terms involving A\\ and A±, it is not clear how to go back from 
to (p in this case. 

Remark. In the proof of Proposition 12. 4^ the quantity A(/? appears and 
one might expect the condition = to be related to limiting Carleman 
weights for Dirac. This can also be seen by writing P^^^p = e'^^'^ Po{hD)e~^^^ 
as A + iB where A = Po^hD) and B = PoCVip), and by noting that 

\\Po,vuf = {{A - iB){A + iB)u\u) 

for test functions u, where {A — iB){A + iB) has full symbol 

(a — ib)(a + ib) H {a — ib,a + ib}. 

2i 
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Here {c, d} is the matrix symbol whose (j, k)th element is X^ILil^jZ) dn^}. A 
computation shows that ^{a — ib,a + ib} = h{A(f)Ii\f, and thus the symbol 
of {A — iB)(A + iB) is nonnegative definite if Aip = 0. It is not clear to 
us if one can exploit harmonicity in proving Carleman estimates for Dirac 
in R". There are interesting recent results related to this approach and the 
semiclassical Fefferman-Phong inequality for systems in the work [26j . 

3. CGO SOLUTIONS, SMOOTH CASE 

In this section we wish to construct complex geometrical optics solutions 
to Cyu = 0, where Cy is given in ([2]), ([3]). It will be convenient to consider 
4x4 matrix solutions (that is, every column of the matrix is a solution) of 
the form 

U = e-p/'^iCo + hCi + ... + h^-^CN-i + h^~^RN). (12) 

This is a WKB ansatz for the solution, where p = if + iip is a complex 
phase function with a Carleman weight, Cj are matrices which correspond 
to amplitudes, and at is a correction term. These are called complex 
geometrical optics solutions because the phase is complex. 

We will give a construction for smooth coefficients (Proposition 13. 1|) , in 
which case there is an arbitrarily long asymptotic expansion of the form (|12p 
where the successive terms have increasing decay in h. This can be achieved 
a (f is a limiting Carleman weight for the Laplacian. In the next section 
we consider Lipschitz coefficients (Proposition 14. 2| ). in which case is the 
linear weight for simplicity. Then the successive terms in (I12p will have only 
a limited decay in h, but we compute sufficiently many terms to be able to 
prove a uniqueness result for the inverse problem. 

Suppose that A,q± G C°°{Q). Writing Pp = eP^^{Po{hD) + hV)e-p/^ = 
iPo{Vp) + hCv, inserting the ansatz (fT2]) in the equation CyU = 0, and 
collecting like powers of h, results in the equations 

iPo{Vp)Co = 0, 
iPo(Vp)Ci = -CvCo, 

iPo{Vp)CN~i = -CvCn-2, 

PpRN = —hCyCN-l- 

We now give a procedure for solving these equations. The first equation 
implies that the kernel of Po^Vp) should be nontrivial. The same applies to 
the kernel of PoiVpf = {"^pfh, which gives the condition 

{Vpf = 0. (13) 

This is an eikonal equation for the complex phase. If ip is given, the equations 
for ifj become 

|V^|2 = |V(^|2, V(p-VV' = 0. (14) 

Assume for the moment that (jl4p is solvable. From (|13p we get ker Pq{SJ p) = 
imPo(V/9), since the image is contained in the kernel and rank(Po(V/9)) = 2. 
We choose Co = PoiV p)Co where Co is a smooth matrix to be determined. 
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Moving on to the second equation, we use the commutator identity 
CvPo{yp) = Ma + Po{Vp){-Po{D + A) + Qj), 
where Qi = ^^j^ ^ and Ma is the transport operator 

Ma = (2(Vp • {D + A)) + -Ap)h. 
The equation for Ci then reads 

Po(Vp)Ci = Po{Vp)-iPo{D + A)- Qi)Co + iMaCq. 

A solution is given by C\ = j{Po{D + A) - Qi)Cq + Pq{V p)Ci, for some Ci 
to be determined, provided that 

MaCo = 0. (15) 

Under certain conditions which are stated below, the transport equation 
(jl5p has a smooth solution Cq. 

For the third equation we use that 

{Po{D + A) + Q){Po{D + A)- Qi) = Ha,w, (16) 

where Ha,w = {D + A)'^!^ + 1^ is the magnetic Schrodinger operator, with 

W = ( ^ • (V X ^) - Q+Q-h -(7 ■ Dq+ \ 

V • ^1- a-(V xA)- q+q-h ) ' 

Then the equation for C2 becomes 

Po(Vp)C2 = Ha,wCo + Po{Vp)-{Po{D + A)- Qj)Ci + IMaCl 

Choosing Ci as a solution of the transport equation MaCi = iHA,wCo, we 
obtain C2 as 

C2 = -{Po{D + A)- Qi)Ci + Po(Vp)C2. 

Continuing in this way we obtain smooth matrices C3, . . . , Ctv-i- The 
equation for /i^~^i?7v may be solved by Proposition 12.31 which ends the 
construction of solutions. 

We still need to consider the solvability of the eikonal equation (jl4p and 
transport equation ([15]). As discussed in [T7], these equations can be solved 
provided that is a limiting Carleman weight for the Laplacian, under 
some geometric assumptions. Suppose that Q CC 0, and that (p is a limiting 
Carleman weight in 0,. Assume that 

is contained in the union of integral curves of V(/J, all 
passing through the smooth hypersurface G = ip~^{Co) in 

n. 

Then one may solve the eikonal equation iVV'oP = I^V^P on G by letting ■0o 
be the distance in the metric iV^Jpeo on G to a point or hypersurface (eo is 
the induced Euclidean metric), chosen so that tpQ is smooth. The limiting 
Carleman condition implies that i/j, obtained from ipQ by extending it as 
constant along integral curves of V(/5, will solve the eikonal equation ()14p . 
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Given ^p satisfying ([H]), one has [Vip, V^p] = cVip + dVip by [T7| (see also 
[TBI Lemma 3.1]). It follows by [8j that transport equations such as ([15]) are 
solvable, provided that the leaves of the foliation generated by Vip and Vip 
are not contained in fi. 

If Lf is the linear weight or the logarithmic weight log|j; — xo|, with xq 
outside the convex hull of 0,, then these geometric conditions are satisfied 
and the equations ()14p and (llSp may be solved explicitly (for the logarithmic 
weight see [?])• Thus, the following result holds in particular for these choices 
of the weight cp. 

Proposition 3.1. Let fhe a limiting Carleman weight for the Laplacian in 
where CC and where the above conditions for solving (jl4p and (jlSp 
are satisfied. If A, q± £ C°°{Q), then the equation CyU = has a solution 
of the form (I12p . where 



Ci = }{Po{D + A)- Qj)Co + Po{Vp)Ci, MaCi = iHA,wCo, 

Cn-1 = \{Po{D + A)- Qi)Cn-2 + Po{Vp)Cn^i, Cn-i smooth, 
llCjlltyi.-(n) <1, \\RN\\+h\\VRN\\<l. 



Here we construct solutions in the case of compactly supported Lipschitz 
continuous coefficients, A,q± G wl'°°{Vl). The construction will be carried 
out for the linear phase p{x) = C, ■ x where C, G C^, = 0, and |ReC| = 



To deal with nonsmooth coefficients, we introduce mollifiers 77e(x) = 
e~^r]{x/e) where 77 G C^(R^) is supported in the unit ball and < 77 < 1, 
Jr]dx = 1. Decompose A = + A^ and Q = + Q\ where = A*r]e 
and = Q *rj^, with the special choice 



where it > is small. Also write V = V'^ + V^ where V'^ = Po{A^) + Q*. We 
find a solution to CyU = of the form 



Co = Po(Vp)Co, MaCo = 0, 



4. CGO SOLUTIONS, Lipschitz case 



ImCI = 1. 



e = h'' 



U 



e-p'^{Co + hCi + R) 



where 



^Po(C)Co 
PpR 



0, 




We make the choices 



C7o = Po(C)e''^", 



4, 
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where cj)^ is the solution to Q ■ (Vcf)^ + ^4") = given by 

<A« = (c-v)-H-c-^«). 

Here {C ■ V)~^ is the Cauchy transform 

(C • V)-V(x) = / -^—f{x-yiReC-y2lraC)dyidy2. 
2vr yi + iy2 

The following result will be used to solve for the error term R. We will need 
the additional small parameter h to prove some extra decay for R. 

Proposition 4.1. Let C R" be a bounded open set, let Pq be a Dirac 
operator in R", and assume that V € W^^"- n L~(17)^^^. Also let C G C", 
= 0. If h, h are small, then for any / G H^{Q)^ there is a solution 
u G of the equation 

(Po(/i-D + C) + ^^)?^ = / inl^, 

which satisfies 

MINI +/^l|Vn||)< 11/11 +/i||V/||. 

Proof. We wish to show the Carleman estimate 

h\\{hD)-'^u\\ < C\\{hD)-^{Po{hD + () + hV)u\\. (17) 

The result follows from this in a standard way by the Hahn-Banach theorem. 

Let p{x) = C-x and pe = P+j^- Let u G C^i^), and choose x G C^{^) 
where x = 1 near CC 0,. Then the pseudolocal estimate 

11(1 - x){hD)-\\\ < CMh^'\\{hD)-\\\, V G ^7^(5^), any M, 

and the Carleman estimate (jlOp imply that when h, h are small, 

h\\{hD)-\\\ < h\\x{hD)'\\\ + h\\{l - x){hD)-\\\ 

< CV^\\Po{hD + Vpe){x{hD)-\)\\ + ^\\{hD)-\\\ 

< CV^WxPoihD + Vpe)i{hD)-\)\\ + ^\\{hD)-\\\. 
Since V/o^ = (1 + " ^)C^ we obtain 

< CV^\\x{hD)-'PoihD + V/9,)^x|| + — ||xPo(C)[C • x, 



Here [C • x, (hD) ^] = hR{hD) ^ where R is bounded on with norm < 1. 
If h is small enough, we obtain 

h\\{hD)-'^u\\ < C^/e\\{hD)-^Po{hD + V/9£)n||. 

If e is small enough, the estimate remains true with Po^hD + Vpe) replaced 
by Po{hD + Vpe) + /iT^, since 

\\{hD)-^Vu\\ < C\\{hD)-^u\\. 

Then (I17p follows since e'''^/'' = me^^^ with ||?7i||y[/i,oo(f2) bounded. □ 



CARLEMAN ESTIMATES FOR DIRAC 14 

Proposition 14.11 and the fact that V G W^''^, gives a correction term R 
satisfying 

\\R\\ + h\NR\\ < \\h/:lCi + v\Co + hCi)\\ + h\\V{hdyCi + v\CQ + hCi))\\ 
< h^-^ + \\V^\\l2 + h^'^'^h + h\\VV^\\L2 + HV^Wl^. 
Choosing h = h"'^ for o"i > smah, one has 

||i?||=o(l), \\VR\\=o{l) 
as /i — > 0. Thus we obtain a solution 

U = e-^(Po(C)e*^" + -{Po{D + A^) - Q\)e'^' h + R), 

where ||i?|| = o(l). Also, since Pp = iPo{Vp) + hCy, we have 

\\Po{C)R\\<h\\R\\Hi+h^\\Cv»Ci\\ + h\\Co + hCi\\L^\\V^\\=o{h). 

Noting that by (llSp we may replace e**^' by where A; • C = 0, we have 

arrived at the solutions for Lipschitz coefficients. 

Proposition 4.2. Let A, q± £ Wc'°°i^), and let C G satisfy = and 
|ReC| = |ImC| = 1. There exists a solution U e to CyU = in 

Q, of the form 

U = e-K-^e*^"+^'=-"(Po(C) + -{PoiycP^ + k + A^)- q\) + R), 
where k G with k ■ C = 0, and where ||i?|| = o(l), ||Po(C)-R|| = o{h) as 

5. Uniqueness result 

We will prove Theorem 11.21 The first step is a standard reduction to a 
larger domain. 

Lemma 5.1. Let Q CC 0,' be two bounded open sets in R'^, and let 
^j^Q±,j £ W^'^^ii}') satisfy Ai = A-i and g±,i = g±,2 in If = 

in then = Cy^ in fi' and 

/ Ul{Vr-V2)Uxdx = ^ (18) 

for any solutions \Jj £ H^n')^''^ of Cy^Uj = in 17'. 

Proof. If Cviu' = in 17', then there is G Ti.{Q)'^ such that vCv2''^ = in 
and v±\Qfi = n'_|_|gn- Let v' = v in Q and v' = u' in $7' \ 17. It is easy to see 
that v' G n{^'f and Cy^v' = in 0', showing that Cy^ C in Q'. The 
same argument in the other direction gives Cv\ = Cy^ in Q'. 
Let f/j be as described. Writing {U\V) = j^, V*U dx, we have 

((^1 - V2)Ui\U2) = -{Po{D)Ui\U2) + {Ui\Po{D)U2). 

Since Cy-^ = Cy^, there is C/2 G 7^(17')^^^ with CyJJ2 = in 17' and also 
{Ui - U2)± G HliQ.'f'^. Thus, writing f/i = {Ui - U2) + U2 and integrating 
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by parts, we obtain 

((^1 - V2)Ui\U2) = -{Po{D)U2\U2) + {U2\Po{D)U2) 
= {V2U2\U2) - {U2\V2U2) = 0. 

□ 

Assume the conditions of Theorem 11.21 By a gauge transformation we 
may assume that the normal components of Aj vanish on dO,. Then by 
Theorem 11.31 we know that Ai = A2 and q±^i = q±^2 on 90. Let Q' be a 
baU such that Q CC 51', and extend Aj and q±j as compactly supported 
Lipschitz functions in $7' so that Ai = A2 and q±^i = q±^2 outside fl. 

Lemma 15.11 shows that in the proof of Theorem 11.21 we may assume that 
Q is a ball, the coefficients are in Wc'°°(0), Cv, = in n, and ^ 
holds for solutions in 57. The recovery of the coefficients proceeds similarly 
as in [22], using now the solutions provided by Proposition 14.21 We will 
give the details since one needs to ensure that the estimates for nonsmooth 
solutions are sufficient for this argument. Theorem 1 1 . 2 1 will follow from the 
two propositions below. 

Proposition 5.2. V x Ai = V x in 57. 

Proof. Choose C G with = and |ReC| = |ImC| = 1. Let k e he 
orthogonal to Re^ and lm(. By Proposition 14.21 we may choose solutions 
to Uj = in 57 of the form 

C/2 = e'K-(Po(C)e^'^" +i?2), 
where (j)] = (C • V)~i(-C • A^-) and \\Rj\\ = o(l) as /i ^ 0. Then 

U^ = e-ii^-^{Po{Oe-^^'^+Rl). 
Inserting these in (fT8]l and letting h ^ gives 

I e^'^-^e'^PoiOPoiAi - ^2)Po(C) dx = 0, 

where (j) = {( ■ V)-^{-C ■ (^i - ^2)) and we have used -Po(C)Qj^o(C) = 0. 
Using the commutator identity for Pq one obtains 

J e'''-''e''^{C-{Ai-A2))dx = 0. 

Lemma 6.2 in |27j, which is based on [lOj, implies that the same identity is 
true with e**^ replaced by 1. Consequently 

J e'''-'^{C-{Ai-A2))dx = 0, 

for all k orthogonal to Re^, Im^. This implies the vanishing of the Fourier 
transform of components of V x (Ai — ^2)- D 

Proposition 5.3. g± 1 = g± 2 in 57. 
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Proof. Since V x = V x A2 and Aj G wl'°°{Vl) where is a ball, one 
obtains Ai — A2 = Vp for p £ W'^'°^{Q) where one may choose p\dn = 0. 
Thus, by a gauge transformation we may assume that Ai = A2 = A where 
A e Wc'^^in). Fix k eK^ and take C G with ("^ = and k ■ ( = 0. By 
Proposition 14.21 we take solutions to Cy^ Uj = of the form 

Ui = eK-^e^'=-^e*^"(-Po(C) + f (Po(V<^" + k + A^) - Q\ j) + 

U2 = e-K"-e^*""(Po(C) + ^{PoiV^i + Ai) - qI j) + R2), 

where = (C-V)-i(-C-^") and \\Rj\\ = o(l), ||Po(C)^i|| + ll-Po(C)^2|| = o(/i) 
as h ^ 0. Also, 

U; = e-K-e-<^"(Po(C) - y (Po(V<A« + A») - + i?^), 
with ||i?2^o(C)ll = o{h). The identity ((IH]) implies 

J e^'^-^mC) - y(Po(V0tt + A«) - g» ,) + - Q2) 

■ (-i^o(C) + -(Po(V0tt + A: + - q! ,) + Ri) dx = 0. 

The highest order term in h is Po{C){Qi ~ Q2)Po{C)j and this vanishes. Also 
the terms involving PqCV^" + go away, since for Q = Qi — Q2, 

Po{C)QPo{y<P^ + + Po(V</.« + A»)QPo(C) 
= Q/(Po(C)i'o(V</.« + + Po(V</.« + ^»)Po(C)) 

= 2g7(c- (V0» + ^«)) = 0. 

All the terms involving Rj are o{h). For the term involving R2QR1, this is 
seen by using and writing 

Ri = l={Po{c)Po{C)Ri + PoiOPoiORi) = PoiCMi) + o{h). 

Therefore, dividing the integral identity by h and letting h ^ gives 

J e^'-^iPoiOQiPoik) - Qij) - Q2,iQPoiO] dx = 0. 

Commuting Po(C) to the left gives 

J e^^-'^PoiOiQPoik) - QiQij + Q2Q2,i) dx = 0. 

This is true also when ^ is replaced by ^, and adding the two identities and 
multiplying by Pq (Re C) on the left implies that 

J e^'^-^HQi - Q2)Po{k) - q+,iq-,ih + q+,2q-,2h) dx = 0. 

Looking at the off-diagonal 2x2 blocks shows 

J e''-''{q±,i-q±,2){a-k)dx = 0. 

The claim follows upon multiplying hy a ■ k. □ 
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6. Boundary determination 

In this section we show that q± and the tangential component of A at 
the boundary are uniquely determined by the Cauchy data set Cy. More 
precisely, we have the following. 

Proposition 6.1. Let Aj, q±j be W^''^{Q) coefficients, j = 1,2, and let Q 
have boundary. If Cy^ = Cy-j , then 

(^1 - A2){xo) ■i = Vxo G Vt G T^oidn), 

Q±,i{xo) = q±,2ixo) yxo £ dn. 

Following the idea of [3j we first construct a sequence of solutions which 
concentrate at xq in the limit. We assume without loss of generality that 
xq = and that Q is defined by the function p £ (R^ ; R) in such a way 
that = {x G R^ ; p{x) > 0}, 90 = {x G R^ ; p{x) = 0}, and the outer 
unit normal of dO. at is —63 = — V/?(0). 

Let now r][x) be a smooth function supported in -6(0, 1/2) such that 

/ r}{x' dx' = I. 
For all M > 0, define r]M{x) := ri{M{x' , p{x))). Define 

where t G To(di}) is a unit vector and is chosen so that 

with Lo{-) a modulus of continuity for Vp. We prove the following. 

Proposition 6.2. For M large enough, one can find solutions to 

the equation 

CvU = 

of the form 

U = NPQ{i + iVp)uo + R 

with 

\\R\\l2 < CN-^'^. 

To prove this result, we will use the following two lemmas from [3j. 

Lemma 6.3. Let rj be smooth and supported in i?(0,l/2). If r]M{x) = 
r]{M{x',p{x))) and 7V"^ = M-^uj{M-'^), we have 

lim M^iV / exp{-2Np{x))'nM{x) dx = 1/2 / r]{x',0)dx' 

and 

I / exp{-2Np{x))r]M(.x)dx\ <CM-^N-^. 
Jn 
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Lemma 6.4. Let A : Q ^ C'^ be a continuous vector field and k : ^} ^ C 
be a continuous function. If is not an eigenvalue of the operator 

-A + A-D + k:H^{n)^ H''^{n), 

then there exist solutions to 

{-A + A-D + k)u = 

of the form 

U = Uq + Ui 

with Ml e i^o(^) and ||ui||/fi(f7) < CN''^/'^. 

At this point we recall some matrix identities used already before. If 
Q = ( ) = ( '~o' g+h ) if a, & e C3, we have 

Po(a)Po(fe) + Poib)Po{a) = 2a ■ b, Po{a)Q = QiPoia). (19) 

Proof, (of Proposition 16. 2p Recall from (jl6|) that 

£y(Po(^ + ^) -Q/) = -A + 2^-Z) + W^ (20) 

where is a 4 x 4 matrix function with entries. Choose A such that 
is not an eigenvalue of the scalar operator 

-A + 2A- D + X: H^{n) H-^{n). 

Apply Lemma 16.41 to get a scalar solution u £ H^{Q) to 

{-A + 2A-D + X)u = 

of the form 

■U = Mo + Ml, 

with Ml G Hq{^) and ||mi||hi < CN^^^"^. We seek solutions to CyU = of 
the form U = {Po{D + A) — Q j)u + R. Plugging this ansatz into the equation 
and using ([^n|) we get that R satisfies 

CyR = (A - W)u. 

By Lemmas 16.31 and 16.41 we see that ||m||/^2 < CN~^/'^. 

If one had well-posedness for the boundary value problem for £y, we 
could solve for R uniquely and obtain the estimate ||-R||2,2 < CN^'^^'^. More 
generally, Proposition 12.31 with the choices (p{x) = xi and h = ho where 
ho = ho{A,q) is the upper bound for h, gives a solution R = e~^^/^°R 
satisfying ||^||j^2 < C\\R\\l2 < C\\e''^/''o {X - W)u\\l2 < CN'^I'^. Writing 

U = Po{D)uo + {Po{A) - Qi)uo + {Po{D + A) - Qi)ui + R 

and using the estimates of Lemmas 16.31 and 16.41 we have the desired form 
for the solution U. □ 

Proof, (of Proposition 16. ip We may assume, after gauge transformations 
if necessary, that the normal components of Aj at the boundary are null 
for j = 1,2. Let Uj (j = 1,2) be solutions to Cy^Uj = constructed in 
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Proposition 16.21 By the assumption that Cy^ = Cyj, we have as in Lemma 
15.11 the orthogonahty condition 

= [ U2iVi-V2)Uidx 
Jn 

= N^ [ Po{i-iVp){Vi-V2)Po{i + iVp)rilie-^^Pdx + 0{M-^). 
Jq 

Multiplying by M'^N~^, taking M — > oo, and using Lemma 16.31 we get that 
at the origin 

= Po{i - iVp)Po{A)Po{i + iVp) + Po{i - iVp)QPo{i + iV p) (21) 

where A := (^i - A2){0) and Q := {Qi - Q2){0). Writing C = i + iVp{0) 
and applying the matrix identities ()19p in ()2ip . we get that 

= -Po{A)Po{OPo{0 + 2{A ■ t)Po(C) + QiPoiOPoiO- (22) 
This is true also when t is replaced by —t, so we have 

= -Po{A)Po{C)Po{C) + 2{A ■ t)Po(C) + QiPoiOPoiO- (23) 

Adding and (^H]) together and using |Cp = 2 we see that 

= -Po{A) + {A-i)Po{i) + Qj. (24) 

If A were nonzero one could take t = ^4/1^41 and obtain Qj = 0, and then 
choosing t G To((9il) orthogonal to A would give ^ = 0. Therefore, ^ = 0, 
and going back to (j24p gives Q = 0. □ 

References 

[1] A. Boutet de Monvel-Berthier, An optimal Carleman-type inequality for the Dirac 
operator, Stochastic processes and their apphcations in mathematics and physics 
(Bielefeld, 1985), Math. AppL, voL 61, Kluwer Acad. PubL, 1990, pp. 71-94. 

[2] R. M. Brown, Recovering conductivity at the boundary from the Dirichlet to Neumann 
map: a pointwise result, J. Inverse Ill-Posed Probl. 9 (2001), 567-574. 

[3] R. M. Brown and M. Salo, Identifiability at the boundary for first-order terms, Appl. 
Anal. 85 (2006), no. 6-7, 735-749. 

[4] A. P. Calderon, On an inverse boundary value problem. Seminar on Numerical Anal- 
ysis and its Applications to Continuum Physics, Soc. Brasileira de Matematica, Ri'o 
de Janeiro, 1980. 

[5] M. Dimassi and J. Sjostrand, Spectral asymptotics in the semi-classical limit, London 
Mathematical Society Lecture Note Series 268, Cambridge University Press, 1999. 

[6] D. Dos Santos Ferreira, C. E. Kenig, M. Salo, and G. Uhlmann, Limiting Carleman 
weights and anisotropic inverse problems, preprint (2008), arXiv;0803.3508 

[7] D. Dos Santos Ferreira, C. E. Kenig, J. Sjostrand, and G. Uhlmann, Determining a 
magnetic Schrddinger operator from partial Cauchy data, Comm. Math. Phys. 271 

(2007) , 467-488. 

[8] J. J. Duistermaat and L. Hormander, Fourier integral operators. II, Acta Math. 128 
(1972), 183-269. 

[9] M. EUer, Carleman estimates for some elliptic systems, J. Phys. Conf. Series 124 

(2008) , 012023. 

[10] G. Eskin and J. Ralston, Inverse scattering problem for the Schrddinger equation with 
magnetic potential at a fixed energy. Comm. Math. Phys. 173 (1995), 199-224. 

[11] , On the inverse boundary value problem for linear isotropic elasticity. Inverse 

Problems 18 (2002), 907-921. 



CARLEMAN ESTIMATES FOR DIRAC 



20 



[12] M. Goto, Inverse scattering problem for Dirac operators with magnetic potentials at 

a fixed energy, Spectral and scattering theory and related topics (Japanese) (Kyoto), 

no. 994, 1997, pp. 1-14. 
[13] H. Heck, G. Uhlmann, and J.-N. Wang, Reconstruction of obstacles immersed in an 

incompressible fluid, Inverse Probl. Imaging 1 (2007), no. 1, 63-76. 
[14] V. Isakov, Carleman estimates and applications to inverse problems, Milan J. Math. 

72 (2004), 249-271. 

[15] H. Isozaki, Inverse scattering theory for Dirac operators, Ann. I. H. P. Physique 

Theorique 66 (1997), 237-270. 
[16] D. Jerison, Carleman inequalities for the Dirac and Laplace operators and unique 

continuation. Adv. in Math. 62 (1986), no. 2, 118 134. 
[17] C. E. Kenig, J. Sjostrand, and G. Uhlmann, The Calderon problem with partial data, 

Ann. of Math. 165 (2007), 567-591. 
[18] K. Knudsen and M. Salo, Determining nonsmooth first order terms from partial 

boundary measurements, Inverse Problems and Imaging 1 (2007), 349-369. 
[19] H. Koch and D. Tataru, Carleman estimates and unique continuation for second order 

elliptic equations with nonsmooth coefficients, Comm. Pure Appl. Math. 54 (2001), 

no. 3, 339-360. 

[20] X. Li, On the inverse problem for the Dirac operator. Inverse Problems 23 (2007), 
919-932. 

[21] N. Mandache, Some remarks concerning unique continuation for the Dirac operator, 

Lett. Math. Phys. 31 (1994), no. 2, 85-92. 
[22] G. Nakamura and T. Tsuchida, Uniqueness for an inverse boundary value problem 

for Dirac operators, Conun. PDE 25 (2000), 1327-1369. 
[23] G. Nakamura and G. Uhlmann, Global uniqueness for an inverse boundary problem 

arising in elasticity, Invent. Math. 118 (1994), 457-474. 
[24] , Erratum: Global uniqueness for an inverse boundary value problem arising 

in elasticity. Invent. Math 152 (2003), 205-207. 
[25] P. Ola and E. Somcrsalo, Electromagnetic inverse problems and generalized Sommer- 

f eld potentials, SIAM J. Appl. Math. 56 (1996), no. 4, 1129-1145. 
[26] A. Parmeggiani, A remark on the semiclassical Fefferman- Phong inequality for certain 

systems of PDEs, preprint (2008). 
[27] M. Salo, Semiclassical pseudodifjerential calculus and the reconstruction of a magnetic 

field, Comm. PDE 31 (2006), 1639-1666. 
[28] M. Salo and ,J.-N. Wang, Complex spherical waves and inverse problems in unbounded 

domains. Inverse Problems 22 (2006), no. 6, 2299 2309. 
[29] J. Sylvester and G. Uhlmann, A global uniqueness theorem for an inverse boundary 

value problem, Ann. of Math. 125 (1987), 153-169. 
[30] G. Uhlmann, Inverse boundary value problems for partial differential equations, Pro- 
ceedings of the International Congress of Mathematicians (Berlin), Doc. Math., vol. 

Ill, 1998, pp. 77-86. 

[31] , Commentary on Alberto P. Calderon's paper: On an Inverse Boundary Value 

Problem, Selecta (to appear) (A. Bellow, C. E. Kenig, and P. Malliavin, eds.). 

[32] G. Uhlmann and J.-N. Wang, Reconstructing discontinuities using complex geometri- 
cal optics solutions, SIAM J. Appl. Math. 68 (2008), no. 4, 1026-1044. 

[33] , Complex spherical waves for the elasticity system and probing of inclusions, 

SIAM J. Math. Anal. 38 (2007), no. 6, 1967-1980. 

Department of Mathematics and Statistics, University of Helsinki 
E-mail address: mikko.salo@helsinki.fi 

Department of Mathematics, Stanford University 
E-mail address: leo.tzou@gmail.com 



